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I. INTRODUCTION 



The effective field theory used most extensively to study QCD at low energies is generically referred to as Chiral 
Perturbation Theory (ChPT). With the inclusion of baryons, the effective theory is called Baryon ChPT (BChPT), 
whose nonrelativistic limit (with respect to the baryons) is referred to as Heavy BChPT (HBChPT). So far, pion- 
nucleus scattering and production processes involving multiple nucleons with an arbitrary number of pions, have been 
considered within HBChPT up to the tree level, with one-loop corrections only at the single-nucleon level (vertex 
corrections) JjLpj- In this paper, we perform a 2-nucleon-l-loop calculation involving pion loops, with the pions being 
emitted and absorbed at different nucleons, which we believe has not been done before. 

The goal of this paper is to determine the size of the 1-loop contributions to pion Double Charge Exchange (DCX) 
scattering at threshold on a nuclear target, relative to the tree graphs, in the framework of HBChPT. One of the 
motivations for this study is the fact that sizable 1-loop (pion) contributions to ir — ir scattering and 7T-N scattering 
have been obtained in the framework of (HB)ChPT by previous authors ( f|||]). So, it's natural to ask whether similar 
large contributions are found for a 2-nucleon calculation involving pions in HBChPT to one loop. We shall find that 
this is not the case; the loop correction to the two nucleon process is small, as expected for a chiral expansion. 

The 2-nucleon process considered is pion DCX: 

7r + + n + n — > 7r~ + p + p, (1) 

where the nucleons are in bound nuclear states. We consider only transitions to the DIAS (= Double Isobaric Analog 
States), e.g. 14 C(7r + , 7t~) 14 (DIAS). The DIAS is that (normalized) state obtained by operating with 7+7+ on 
the target ground state, where 7 is the total isospin operator. The DCX contribution is much smaller than elastic 
scattering ( 14 C^ 14 C) because while the former involves only the valence nucleons, the latter involves a coherent 
scattering of the core and the valence nucleons. 

The reason for considering DCX at threshold is that it is dominated by two-nucleon processes, since single ir — N 
scattering cannot contribute. It has been shown that the contribution of meson exchange currents (MEC) to DCX, 
although less than double-scattering, is not all that small. (See || and references therein.) In the context of HBChPT, 
these results are at the leading order, or tree level. The second motivation of the present paper is to establish the size 
of the next-order correction to the two-nucleon MEC contribution to DCX. 

In Section II, we introduce the theoretical background of single- and multi-nucleon HBChPT and discuss a chiral 
power counting rule due to Weinberg (including issues of reducibility of graphs, which is an additional complication 
arising at the multi-nucleon level). In Section III, we discuss the approximations involved in getting (analytical and 
numerical) estimates of the amplitudes for the tree and 1-loop graphs for pion DCX. The relevant leading order tree 
graphs are evaluated in the framework of HBChPT and compared to an earlier calculation. In Section IV, we evaluate 
the 2-nucleon 1-loop corrections to the tree graphs of Section III. We then make numerical estimate of the finite parts 
of the (2-nucleon) 1-loop graphs, and make comparison with the tree graphs of Section III. In Section V, we discuss 
the renormalization of the 1-loop graphs of Section IV using 2ir - 2nucleon contact terms. Section VI summarizes 
and discusses the findings for the DCX problem. We include an estimate of the effect of vertex corrections on the 
tree graph amplitudes, based on earlier published work. There are two appendices: Appendix A presents the vertices 
written in terms of n ± ,ir° [rather than their cartesian counterparts (as in as the former are more readily useful 
for calculation purposes] (A.l), and combinatoric factors (A. 2) for the 1-loop 7T-NN graphs. Appendix B discusses 
the various 1-loop integrals and related identities involving them, relevant to the DCX 1-loop calculation, highlighting 
the ones that are new. 



II. THEORETICAL BACKGROUND 



In this section, we discuss the basic elements of HBChPT at the single and multi-nucleon level that will be required 
later in the paper. 

The leading order (LO) HBChPT Lagrangian that will be used in the calculations of the LO tree and and 1-loop 
graphs is given by: 



H(iv ■ D + g^S ■ u)H 



(2) 



where g\ = axial- vector coupling constant, F = pion-decay constant and M =pion mass in the chiral limit. The trace 
in the nucleon isospin space is denoted by ( ) in (||). 
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The HBChPT Lagrangian is written in terms of the "upper component" H and its covariant adjoint fi, exponentially 

parametrized matrix-valued meson fields U, u = VU, baryon ("u^S,,") and pion-field-dependent ("D^, u u , x±") 
building blocks defined below: 

H = e w«l (1 + jW) (3) 

where "0 is the Dirac spinor and m is the nucleon mass; 

= nucleon 4 — velocity parameter, 

S„ = -7 5 Ci/ P w p = Pauli — Lubanski spin operator; (4) 



U = exp i— I , where </> = 7? ■ t, (5) 

V -Fir/ 

where t are the nucleon isospin generators; D M = <9 M + where = \[u\ d^u]; = i(w d^u — ud^u^). 
The run — NN vertices, m — 1, 3, are constructed from the Yukawa term: 

.9° HS • u^H, (6) 

(where the superscript on represents the powers of the pion field). The 2-k — NN vertex is constructed from the 
Dirac term: 

iEv • D (2) H = iRv • T (2) H (7) 

(where the superscripts on D^, and represent the powers of the pion field). 

The four-pion vertex is constructed from the non-linear sigma model Lagrangian = LO ChPT Lagrangian: 

" ' • (8) 
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■ ({WUidnU) + M 2 (U + - 2)\ . 



For more details, refer to Appendix A. 

The elementary vertices from which the tree and the 1-loop graphs of Sections II and III have been constructed, are 
drawn in Fig 1. However, the calculations in the chiral perturbation expansion are renormalized to the chiral order 
of the expansion. The Weinberg chiral power counting relation (WCPCR, is used for a systematic classification 

of the relevant tree and 1-loop graphs. The relation determines the overall chiral order of irreducible graphs in terms 
of the total number of incoming or outgoing nucleons (the two are the same because of baryon number conservation) , 
the total number of loops, the chiral order of the vertices, and connectedness of the graphs, as discussed below. Here 
is WCPCR: 

v = A-N + 2{L-C)+Y J ViC^+d i -2), (9) 

i 

where v = overall chiral order of a graph, TV = total number of incoming/outgoing nucleons, L = number of loops, 
C = number of separately connected pieces of the graph, Uj = the number of vertices of type i, rij = the number 
of incoming and outgoing nucleons at the ith vertex and di = is the number of derivatives or powers of M v . From 
baryon number conservation, m = 2 or 0. 

Graphs which violate relation (^) because of anomalously small denominators are referred to as reducible by 
Weinberg pi. This class includes graphs whose energy denominators in old-fashioned time-ordered perturbation 
theory are of the order of M„ /m rather than M v . These graphs arise in the context of NN or many-nucleon 
scattering, without external pions. (For a recent treatment of the NN scattering problem which deals with problems 
of the Weinberg scheme, see pjLQj.) However, for pion scattering on NN pairs, as in the present paper, all graphs 
are irreducible, as can be seen from the discussion of ||. In Sections III and IV, WCPCR will be used to determine 
the overall chiral order of the tree and 1-loop graphs. 

The lowest order tree graphs for DCX, based on the LO vertices @-(||), will be shown to be of chiral order v = 
(see Sec. III.B). The one-loop corrections at the 2N level, which are the main subject of this paper, are of order v = 2 
(see Sec. IV). To this same order, there are also corrections to the LO vertices, which have already been obtained in 
(HB)ChPT (see ||[|,|ll]]) as renormalized effective interactions, with a number of low energy constants (LECs). These 



3 



LECs are to be fixed from experiment, but are not in fact all known. (For the tt — tt vertex, the LECs of v = 2 are 
'almost' all determined [9. For tt — N, the information is less complete, but has been supplemented by theoretical 
arguments for on-shcll nucleons.) 

However, we do know the renormalized values of the pion and nucleon masses, the axial- vector coupling constant, 
and the pion-decay constant. For the purpose of numerical calculation of the analytical expressions for the tree and 
1-loop graphs later in the paper, we use 

A'U = 139.57 MeV, F n = 93 MeV, g A = 1.26. (10) 

Renormalizing these constants, including the nucleon mass to, to 1-loop (in the tt — N and tt — tt interactions) gives 
the following relations to the corresponding quantities in the chiral limit g A , to , M, F : 

g A = g° A [l + Pg M 2 } 
m = m°[l + p m M 2 } 
M 2 = M 2 [l + p M M 2 } 

F n =F[l + p F M% (11) 

where the pi include some of the (HB)ChPT LECs. 

It is also known that the 1-loop corrections to the tt — tt vertex give contributions of about 25% of order v = 2. 
Similar corrections to the to,7t — N vertices (with m — 1, 2, 3) are as large as 10 - 30% ||Q,[ll]]. So we know the order 
of magnitude of the v = 2 corrections to the LO DCX (tree) amplitudes from the 1-loop vertex contributions, but 
cannot determine these corrections completely, without the undetermined LECs. 



For that reason, we shall use the vertices (|6|) - (g), with the renormalized values (10) for the tree calculation in 
Sec. Ill, omitting unknown (but significant) corrections of the same order. The 1-loop corrections at the 2N level are 
calculated in Sec IV with the same vertices. 

III. DCX SCATTERING AMPLITUDES; TREE GRAPHS 

In this section we discuss the approximations that will be made in evaluation of tree and 1-loop graphs for DCX 
scattering of pions off a nuclear target, set up the notations and calculate the amplitudes for the leading order (DCX) 
tree graphs. 

A. Notations and Approximations 

The scattering matrix element Sfi is defined as: 

S fi = -i(2n)W\P f -Pj-^-L—M. (12) 
In (|l2]), the nuclear scattering amplitude M. is defined as the matrix element of a two-body operator T: 

M = (iJj , 1=1, h = l\Tr^T^\tP h 1=1, h = -1) 
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= <Vo|T|^) x (/=l,7 3 = l|r| 1) r| 2) |7=1^3 = -l), (13) 



where we assume a target with 2n+ isoscalar core (I = 1, 13 = —1); then the DIAS has 2p+ isoscalar core (I = 1, 13 = 
1). The isospin matrix element in the second line of (^) equals unity, and will be omitted in the following. 

If 2 were the 4- momenta of the incoming nucleons, pt^ 4 the 4- momenta of the outgoing nucleons, and q^ 2 the 
4-momenta of the incoming and outgoing pions (respectively) then 

r 4 

M = J Y[d 3 Pl (MP3,P±)\T\A(Pi,P2)). (14) 

i— 1 

The nucleons are in a relative I = state (1 = state being the dominant partial wave for ground states) and hence 
from Pauli's exclusion principle, the spin of the incoming and outgoing nuclear states must be equal to zero. Hence, 
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from here on, for both tree and 1-loop graphs, we will simplify the structure of the transition operator T assuming 
that eventually one is going to take it's expectation value with respect to \l = S = 0) nuclear states. 
The following will be used extensively in the same context. 

(S = 0\(a^ . $(SW . g)\ S = 0) = • a (2) )s=oy = (15) 

Following (lj and First, the velocity parameters of the two participating nucleons are both chosen to have the 
static limit values, with only a non- vanishing time component , i.e. v 1 ^ = = (1,0). Second, the nucleons will be 
treated as if they were on-shell (sometimes referred to as impulse approximation). In the HBChPT formalism, p° 
denotes only the contribution to the time component of the total nucleon momenta (= mv + p) , in addition to the 
rest mass energy m (for the choice of the nucleon velocity to possess only a non-zero time component). In the present 
case p° = Eb, which as stated above, we drop. Then if we go to the cm. frame of the nucleons: 

= (0,p); p» = (0,-p); p% = (0,if); pi = (0,-^); <# = g£ = (M^O). (16) 

(Note: The external pions are at zero kinetic energy [threshold].) 

For this paper, we shall only evaluate (S = 0\T\S = 0). M will be estimated by calculating (S — 0\T\S = 0) at a 
typical point: 

P 2 = M%, (17) 

where P = p 1 — p. This is a reasonable kinematic point because the inter-nucleon separation in a nucleus averaged 
over the nuclear wavefunction is roughly M~ l . Then 

r 4 

M « (s = o\t(p 2 = Ml)\s = 0)1 Y[d 3 Pl (MP3; PiMtipi; P2)) 

= (S = 0|T(P 2 = M 2 )\S = 0) (18) 
where the last line follows because the overlap integral is unity for ip = DIAS of ipi. 



B. Tree Graphs 

The connected tree graphs are of overall O(q ), as follows from the chiral counting law (^|), because for them 
iV = 2, L = 0, G = 1, Ui = 2 or 0, di = 1 or 2. The tree graphs have been constructed from the first three and the 
fifth elementary vertices of Fig 1. They are drawn in Fig 2. For LO graphs, all the vertices have to be of LO, i.e. 
O(g) mi: — NN (m = 1,2,3) and 0(q 2 ) 4 ir vertices. Using the exponential parameterization for the matrix-valued 
meson field, /ChbciiPT and £chPT are written down explicitly in terms of the pion triplet fields : tt°, ^ and p, n fields 
[refer to Appendix A.l]. 

The S-matrix for (Q) can be calculated in perturbation theory using standard Feynman-diagram techniques. For 
the nth order term , , a combinatoric factor / is defined via: 

' J fc=l H=l j ' 

/n ,n \ 

l[d 4 x k T [HjCfixi) , (19) 
fc=l H=l ' 

(T = time-ordering operator), where one uses: T(AB) = T(BA), where A,B= bosonic fields or fermionic bilinears. 
One can show that the combinatoric factors for tree graphs of Fig 2 are: 

contact graph Fig 2a : f = 1; 
Pole graph Fig 2b f = - ; 

double — scattering graph Fig 2c : f = — . (20) 
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The amplitudes for the tree graphs are written in terms of Pauli spinors (to which H and H reduce in the static 
limit) and the pionic held <j>~J+ (See (^4|)). 

Using also (|l6|), ( |l7| ) and (J2Q), the amplitudes for the tree graphs are expressed as operators written in terms of 1 
and <?« • Pd^> ■ P. 

What follows are expressions for T and (S = 0\T(P 2 = M 2 )\S = 0). 

(a) Contact graph Fig 2a: 



T [a] =2(WJ*-/ 1) - P * 2 l- P . (21) 
N v ; 24F 4 M 2 + p2 ' v > 



which using (||) and @, gives: 



■,2 



(5 = 0|T [q] (P 2 = M^)|5 = 0) = = -1.42 x l(T 8 MeV- 4 . (22) 

(b) Pole graph Fig 2b: 

r^^'^pM^-p r'^y . (23) 

which using @ and @, gives: 

(S|T [fc] (P 2 = M *)|S = 0) = = -6.3 x 10- 9 McV- 4 . (24) 

(c) Double-scattering graph Fig 2c: 



which using (17), gives: 

(S = 0\T [c] (P 2 = M 2 )\S = 0) = = -2.67 x l(T 8 MeV- 4 . (26) 

Thus, the total tree-graph amplitude is : 

(S = 0\T total [tree}(P 2 = M 2 )\S = 0) = - ^ + 2^ = -4.81 x lO^MeV- 4 . (27) 

An older form of chiral Lagrangian that predates QCD was given by Olsson and Turner (OT) ( Q), and has 
been used for tree calculations, e.g., of DCX ||. It consists of the minimum number of derivatives of the pion 
and nucleon fields, with two undetermined (model dependent) parameters. It can be easily shown that the OT 
effective Lagrangian is equivalent to the LO vertices obtained from (jfj) and (^), with the OT parameters taken to 
be J7"=(|, — i). Therefore, the tree graph calculation of this section should agree (as they indeed do) with the 
expressions obtained in || for the transition operators for "forward scattering" for the "pion-contact" and "pion- 
pole" graphs, using the same values of the OT parameters. This provides a check on the tree-level calculation using 
HBChPT. 



IV. 1-LOOP GRAPHS FOR tt-NN 



In this section, we discuss how to evaluate one-loop corrections to the tree graphs (that were evaluated in Section 
III), for pion DCX, in the framework of HBChPT. Loop graphs, in the present context, involve the emission of 
two or more pions from one nucleon leg and their absorption at another nucleon leg (for multi- nucleon processes). 
For evaluation of the 1-loop integrals that occur in the 2-nucleon-l-loop graphs, use has been made of dimensional 
regularization in which the space-time dimension "d" is allowed to vary continuously, and expressions obtained after 
integration, are expanded around d = 4. 
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One gets eight 1-loop graphs using all five of the elementary mir — NN (m = 1,2,3,4) and 4-7r vertices of Fig 1. 
They are drawn in Figs 3 and 4. The uncrossed counterpart of Fig 3d, is not included because it belongs to the class 
of multiple-scattering graphs involving at least three nucleons, which are not considered in this paper. For Fig 4h, if 
the two pions exchanged were 7r°'s, then the amplitudes would vanish in the static limit. 

The leading order connected 2-nucleon 1-loop graphs (for DCX) are of 0(g 2 ), as can be seen from (g), because for 
these graphs, N = 2, L = C = 1, (rii, di) = (2, 1) or (0, 2). The vertices corresponding to the interaction of an even 
number of pions (two or four for our purpose) with nucleons, are obtained from the Dirac term (Q), while the vertices 
corresponding to the interaction of an odd number of pions (one or three for our purpose) with nucleons, are obtained 
from the Yukawa term (g). The LO(=0(g 2 )) ChPT La grangian rtq), which is equivalent to the non- linear a model, 
is used for the 4tt vertex. However, we do use the renormalized constants (|l0|), which introduces some corrections of 
higher order. For details refer to Appendix A.l. 

The one-loop amplitudes are written using standard Feynman rules. The expressions for the T-matrix elements for 
the eight one- loop graphs for fixed momenta of the external nucleons legs, are given below. For notational convenience, 
Hnsntron fai .2) is represented as n(pi. 2 ), and H prot on(P3,4) is represented as p(p3, 4 ). Including the combinatoric factors 
"/" of (|l9|), one arrives at the expressions below for TL-i, j — a,...,h. They are first written in a covariant notation 
except that the two velocities of the two nucleons are chosen to be the same, i.e. v\ — v 2 = v (as in Section III). 
Then the static limit kinematics (|l^) is applied. (The combinatoric factors are given first for each diagram, and are 
obtained in Appendix A. 2). 

(a) 



(V2fg 2 A 1 t d d k 
962* % J {2n) d 

p(p 3 )v ■ (2p 2 - 2p 4 )n(p 1 )p(p 4 )S( 2 ) 



(p4-P2 + k)S^ -fcn(p 2 ) 



(v ■ (fc-pa) - ie)(M 2 - k 2 - ie)(M% - (pi-p 2 + k) 2 - ie) 



(28) 



Using (M 
(b) 



vanishes because v ■ (j>2 — Pa) — * 0. 



d d k 



96F6 i J (2ir) d 

p(p 3 )S( 1 ) ■ fcS^ 1 ) ■ (gi - g 2 + 2 Pi - 2p 2 + 2k)n{ Pl )p(p 4 )v ■ (p 4 - P2 + 2fc)n(p 2 ) 
(v ■ (k- ps) - ie)(M 2 - fc 2 - ze)[M 2 - (k +p 4 - p 2 ) 2 - ie] 

pfe)S (1) ■ (gi -g 2 + 2p A - 2p 2 - 2fc)S< 1 ) ■ kn( Pl )p(p 4 )v ■ (p 4 - p 2 - 2k)n{p 2 ) 
(y ■ (k — pi) - ie)(M 2 - k 2 - ie)[M 2 - (k + p 2 - p 4 ) 2 - ie] 



(29) 



(c) 



(V2r 9 2 A M 2 1 



d d k 



96F6 i J (2n) d 
P(P3)S {1) ■ (Ps ~Pi + k)S^ ■ knipx) , p( P3 )S( 1 ) ■ fcS(D • (pa - Pi + k)n( Pl ) 



(v ■ (k- pi) - ie) 



(v ■ (k+p 3 ) - ie) 



[4gi ■ q 2 + 2k ■ (gi - q 2 + k)]p(p 4 )v ■ (gi - g 2 + P3 - Pi + 2fc)n(p 2 ) 
[M 2 - (k + qi- q 2 ) - ie](M 2 - k 2 - ze)[M 2 - (k - [ Pl - p 3 ]) 2 - ie] 

(d) One can show that the contribution of (d)(1) and (d)(2) (in Fig 3d) are equal, giving an overall factor of 2 

2 5 2 ( % /2) 6 l f d d k 



(30) 



T] 



[d\ 



3 256F6 i 
P{P3)v ■ (qi 



(27r) d 
P2 + Pi 



k)v ■ (k - g 2 )n(pi)p(p 4 )u • (p 2 - p 4 + 2fc)n(p 2 ) 



(v(k + q 2 - pi) - ie)(M 2 - fc 2 - ie)\Ml - (k - [p 4 - p 2 ]) 2 - ie] 



(31) 
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(e) 



(f) 



(g) 



00 



(V2) 8 1 f d d k 



384FJ? 



(2tt)° 



ji ■ q 2 + 2k ■ (qi - q 2 + k)] 



{Ml - fc 2 - ie)[M* - {qi - q 2 + fc) 2 - ie] [M 2 - (k - fa - p 3 ]) 2 - 
P(P3)« ' (pi -P3- 2k)n(p 1 )p(p 4 )v • (qi - q 2 + p 3 ~ Pi + 2k)n(p 2 ) 



If] 



, (V2) 6 1 



d d fc 
(27? 



pw • (g 2 - qi + 2p 2 - 2pi + 2fc)n(pi)p(p 4 )u • (p 2 - p± + 2fc)n(p 2 ) 



{Ml — k 2 — ie)[Ml -{k- [p 4 - P2]) 2 - ie] 



-bl 



g 2 A (V2) 6 1 f d d k ( p(p 3 )S« ■ (pj - p 3 + < Zl )n(p 1 )p(p 4 )S( 2 ) • (pj - p 3 + gi)n(p 2 ) 



144f 6 i 
p(p 3 )S« • (p; 



(27r d ^ (Af 2 - fc 2 - ze)[M 2 - (fc - [gj + Pl - p 3 ]) 2 

pi + 3fc)n(pi)p(p 4 )S( 2 > • (2p 2 - 2p 4 - 3q 2 + 3fc)n(p 2 )\ 



(Ml 



ie)[M%-(k-[ qi + Pl -p 3 ]) 2 -ie\ 



[h] 



, gj(V2) 6 1 r d d k 

'p(p 3 )S« • 



(27r) d 
- P4 + Ql 



92 



3fc)n(pi)p(p 4 )S (2) • kv ■ (fc + p 2 -P4 - 2q 2 )n(p 2 ) 



^ (w • (fc - p 4 ) - ie)(Ml - (fc - [p 4 -P2 + 9 2 ]) 2 - «e](M 2 - fc 2 - ie) 

p(p 3 )S (1) • (p 2 -p± + q\-q 2 - 3fc)n(p 1 )j5(p 4 )t> • (p 4 - p 2 - 2ft + fc)S • kn(p 2 ) 



-<7i 



(«•(*- P2) - ie)(Ml - fc 2 - ie)[Af 2 - (fc - [ 9l + p 2 - p 4 }) 2 

~q2 



ie 



(32) 



(33) 



(34) 



(35) 



The amplitudes for the eight 1-loop graphs are now written in terms of the 11 of the 13 1-loop integrals discussed 
in equations (J51) — (i?3) of Appendix B. They, like the tree graphs of Section III, are evaluated at (|l7|). The Pauli 
matrices o^ 1 )^ 2 ) come from the non-relativistic reduction of S^'' 2 '. The expressions below are written in terms of 
Pauli spin operators, and will be evaluated for (S = 0|T|5 = 0). For graphs (a) - (f), the spin operator is unity; for 
(g) and (h), we use (15). 

Using dimensional regularization, the "L" in the expressions below is the UV-divergent portions of the loop integrals 
and is defined by : 



,d-4 



167H 



l + r'(l) + in(47r) 



(36) 



where d is the dimension parameter and /i is the renormalization point, which is taken to be of the order of the nucleon 
mass. Following jflj we use fx = 1 GeV, and also vary the value by ±0.5 GeV to test the sensitivity of the results to 
this choice. 
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The momentum loop integrals are then rewritten in terms of the 1-loop integrals of (Bl) - (B3): the 2ir - propagator 
integrals are represented by J^,J\2S an( i the 1-nucleon - 2ir - propagator integrals are represented by 70,1,2,3,4,5- 
[] For convenience, they are represented by the 14 integrals U, i = 1 ... 14 in the one- loop amplitudes fivj\) - (|47|). 
Then using the results of Appendix B.l — B.2, the integrals and hence the amplitudes are written down as a linear 
combination of the UV-hnite (L-independent) and UV-divergent (L-dependent) terms. In Table I we identify the 
equation in Appendix B used to evaluate each of the 14 1-loop integrals. The numerical values of the UV-finitc 
parts are also given in Table I, for completeness, so that the calculation can be reconstructed by the reader, or for 
application of the relevant parts to problems other than DCX. 



TABLE I. List of 1-Loop Integrals in ^ - @ 



Ik 


Integrals of App B 


Equation 


Numerical value 
of UV-finite part 
at (i =1 GeV 


h 


J™{-M*,Mi) 


(B9) 


0.018 


h 




(B19) 


0.009 


h 




(B23) 


0.014M 2 


h 




(B25) 


0.006 


h 


QjjrJi^i-MZ, M 2 ) M 2 =M 2 


(B21) 


-0.004A/" 2 


h 


-gjtfj J™ : {—M£, M*)\ M i =M 2 


(B27) 


0.002 


h 


-qj?z Js w (- Mi ,M 2 )\ M 2 =M 2 


(B29) 


-0.004 


h 




(BIO) 


0.019 


h 




(B8,B17) 


0.013 


ho 




(B10,B18) 


0.006 


hi 


7o(-A/ 7r ,0,-M^) 


(B13) 


0.019A/- 1 


I12 


7 2 (0,-M w ,0)-7 2 (0,iW x ,0) 


(B31) 


0.025M- 1 


ha 


73(0, —M-x, 0) — 73(0, M w , 0) 


(B34) 


-0.022M T 


In 


75(0, -M n , 0) -75(0, M«,0) 


(B37) 


0.022M- 1 



x The notations of |^| and [Q, where similar integrals have been evaluated, have been modified for this paper; see Appendix 
B. 
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T [b] = 2(V2) e 



= 2(V2) 6 
29 



48F 6 
9\Ml 



-3/ 3 + Mlh 
29 



1927T 2 



In- 



48F6 
M 2 



576tt 2 + 8^ 2 V\/5-l 



96F6 



-6M< 



-M 2 J 5 - 3/ 6 + Mlh 



+2^Mlh + 3I 3 - Mlh 
g\Ml(V2f 



96FJ 



23 



3 288tt 2 120tt 2 



41-v/5, /\/5 + l 
in 



25 M 2 ^ 
m- 



V5-1/ 48^ 2 m 2 



, (V2) 6 1 
' 3 128^6 

(V2) 6 M 2 
3 64F6 

V5 



2tt 



-in 



1 



-J 3 - IMlh + 4M3j n 



— L 

6 



55 , M 2 53>/5, /V5 + 1 



173 

576tt 2 ' 1927T 2 ""' n 2 ' 1927r 2 """VV5- 1 



-(V2) 



48F6 
M 2 (V2) 8 



h + 3M 2 / 6 



48F6 



3 



1 



96tt 2 



96^6 

(v^fM 2 
48F| 

575 / 
192^ m VV5-l 



1l. 

6 

i + 1 



1 



576tt 2 



M 2 

-29 + 21Zn^L 
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T [g] = 



(V2) 6 gj 



h + 9/io W (1) ' Po {2) ■ P 



(43) 



Using 



one gets: 



One thus gets: 



/io = ^s; / 9 = -^a ot (o,m2), 



~ l9] " 576F6 



(S = 0|T w |S = 0> = -(^. 



1 / 71/ 2 

-24X + ^f3-24/n 



(44) 



(45) 



(46) 



(h) 



[h] 



2 x 



2(V2) 6 gj 
384FJ? 



+ ^-/ 8 - 2M 7r / 12 + 6/10 + 6M X / 14 ) • Pa( 2 ) • P 



2 x 



2(V2) 6 <& 



+ -2L 



384F6 
3 



-6L 



16tt 2 ^ 2 4tt 2 



1 r ^ + A)^(D.p^).p 



16tt 2 16tt 2 ^ 2 32 



(47) 



One thus gets: 



<S = 0|T [fc] |S = 0) = (v^) 6 ^ 



(48) 



Finally, the total 1-loop contribution to the DCX amplitude is given by : 

(S = 0\T total (l - loop)\S = 0) 



24F6 \ 
1 , M 2 



', 2 15 



15 



16 



5 5 j 



181 



+ 



256tt 2 64 
V5 , /V5 + 1 



5a 



24^ 2 VV5-1 



157 
1536tt 2 



11 2 39 
10 5A+ 64 



1 

167T" 



In 



1 + \/5' 



(49) 



The numerical values of the UV-finite parts of the 1-loop amplitudes and their UV-divergcnt portions, are given in 
Table II. 
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TABLE II. Amplitudes for 1-loop graphs 



(S = 0\T [X] \S = 0) 


UV-Finite Integral (units of MeV~ 4 ) 


UV-Divergent Portion 
(coefof^^)L 


a 






b 


-7.7 x 1(T 10 


2 x ^Sl 

288 


c 


1(T 9 


-2 x 112a. 

z 288 


d 


1 x io -11 


sr^ 

1152 


e 


-2.9 x 1CT 1U 


1 

9 


f 


7 x IO -11 


V 

288 


g 


1.9 x ltr 1 " 


9 A 
24 


h 


-10- 9 


5^ 
24 


total 


-(7.9 ±1.6) x 1(T 10 
for fj, = 1 ± 0.5 GeV 


21 x (i+4 fl i) 
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Hence, comparing the numerical values of (|27|) and ( [49D one sees that the 1-loop graphs contribute an increase 
of about a 1.6% relative to the tree graphs, after removal of divergent terms (See Section V). This agrees with the 
expectation that the second-order correction in the chiral expansion should be smaller by a factor of M%/(4ttF w ) 2 ~ 
0.014 than the leading order. The suppression of the 2-nucleon graphs relative to the leading order tree graphs can 
be seen explicitly by comparing the expressions (ETT) and (|49|). (This is in contrast to the large 1-loop corrections for 
the ir — ir and it — N vertices, which is discussed in Section VI.) Also note the relative insensitivity of the total (finite) 
1-loop amplitude to the choice of the renormalization point /i, as shown in the last line of Table II; a 50% change in 
fi gives a 20% change in amplitude. 

From Table II, one sees that the dominant 1-loop contributions come from graphs (b), (c) and (h). This is probably 
related to the following observations. First, as is clear from Table I, the l-nucleon-27r-propagator loop integrals 
dominate over the 27r-propagator integrals. |^ Second, graphs (b), (c) and (h) are proportional to g\ > 1, as compared 
to graphs (d), (e) and (f) (which do not have a g\ dependence). 

V. RENORMALIZATION BY CONTACT TERMS 

In this section we discuss the renormalization of the 1-loop graphs evaluated in Section IV. To remove the UV 
divergence in 3+1 dimensional space, one has to look for terms whose contributions precisely cancel the coefficient of 
the "L" in the (2-nucleon) 1-loop amplitudes. Given the structure of the UV-divergent parts of the 1-loop amplitudes, 
2tt - 2nucleon contact terms of 0(q 2 ) should do the job as the correct counter terms, as will become clear from ( |5l| ) 
and (|52j ) below. 

The contact terms are defined by: 

lim(a'|H0 1 H(x)|a)(/3'|H0 2 H(y)|/3), (50) 

x—>y 

with a = p, I3 . From the expressions for the amplitudes for the eight 1-loop graphs in Section IV, before setting 
P 2 = qi ■ q2 = M 2 , one can show that the UV-divergent terms have the following structures (omitting the overall 
isospin factor t^t^ from (|5ll)) 

/ P 2 or qi ■ q 2 or M% \ x 1; 

P 2 or qi ■ q 2 or M 2 ^j x <r (1) • ct (2) ; 
ffW . pd&) . P. (51) 
On taking the Fourier transform of the amplitudes (|5l|), one gets the local forms: 

~ (Ml, Ql ■ q 2 , dld 2 )S^(xi - x 2 ) = 0(q 2 ). (52) 

Thus one requires O(g 2 ,0 2 ) 2 tt-2 nucleon contact terms with 7+7+ isospin structure. 

For the purpose of renormalization, one needs to consider nine 2-7T - 2nucleon contact terms, written in terms of: 

Up, d p >X+> and 0) ( 53 ) 

Dp = (<9 M + T M ) +(<9 — r^); x+ — M 2 (U + U^)^j. The finite parts of these nine terms will not contribute, as 

explained later in this section. 

One needs to include at least <f> explicitly as a building block for DCX for the following reason. The <f> in (A. 2) can 
be expanded in terms of the generators of the nucleon isospin group as : 

<^ = T30° + V2(t + 4>- +r_0+), (54) 



2 Graph (d) also has a l-nucleon-27r-propagator structure, but the corresponding loop integrals occur in such a way in the 
amplitude that there is large cancellation between the contributions of Ik, k = 1, 3, 11. 
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where r± = i(ri ± 1x2), and </> + either anihilates 7r or creates ir + , and either anihilates ir + or creates it . In 
DCX, a 7r + goes over to 7r~, implying that one requires Oi(x) and 2 (y) in (50) to consist at least of <j>-{x)(j)-{y) (or 

derivatives thereof), or in terms of the isospin generators, t^t+ \ where the superscripts refer to nucleons 1 and 2. 
Further, using Table III: 
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TABLE III. Properties of Building Blocks 



Building Block 


Isospin nature 


Chiral Order 


Dm 


isoscalar(= 9 M ) 
+isovcctor(= r p ) 


O(q) 


u M 


isovcctor 


0(9) 


X+ 


isoscalar 


0(9*) 




isovcctor 


0(1) 



one can construct Table IV for Oi(x)0-2(y) relevant to DCX (again omitting t+V+ 2 ^) which has to include terms of 
the type isovectorxisovector. 



TABLE IV. 0!(x)0 2 (y) relevant to DCX 



Momentum-spin form 


Coordinate-spin operators 


M 2 x 1 
qi ■ qi x 1 

P 2 x 1 


u M (a:)u M (?/); w • u{x)v ■ u(y) 
4>{x) *D,T (x)0(y) (y) 


Ml x a« ■ a< 2 > 
91 • Q2 X • a (2) 

P 2 X • a (2) 


S^V(^)S( 2 )-^(y)(x + } 
« M (a;)SS 1) w' i S I/ ' (2) ; «; • w(a;)S' l ' (1) t; • it(j/)S<, 2) 

0(a;)S (1) ^ *dT (x)0(y)S^ 2) "d" (y) 




#e)S (1) - *~lT (a;)0(y)S (2) - "iT (y) 
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Note that S • u(x)S ■ u(y) will not contribute in the static limit and at threshold (use ( |l6| ) and v± 2 — (1>0)). The 
relevant counter terms are listed below: 



lim 

x—>y 



>(#H(x)H( y )f)H( y )( X+W ) 
+ ||H(x) UM ( a ;)H(x)H(y)u^(y)H(y) 
— jH(x)u • u(x)H(x)H(j/)« • u(j/)H(y) 

^H(,)« DV (,)H(,)H( y) M ^DV (yJHfo) 

+ HSW^^)H(x)H( y )^S( 2 )-H( y )( X+ (x)) 
+ |iH(x)^(a;)SWH(x)H(y)u^(y)S( 2 )^H(y) 
+|^H(a> • uS^R(x)E(y)v ■ u(y)S^U(y) 



§H(,)SW^'D'''(«H(y)S 



+||H(x)«S«. ^ (x)H(x)H(y)M S ( 2 ). ^ (y)H(y) 



(55) 



The coupling constants have been made dimensionless by construction of the terms. 

It should be noted that the following five types of UV- finite 27r-2nucleon 0(q 2 ) contact terms will also contribute 
to off-threshold DCX: 



i^ x v p Ii(x)u^xMx)fl(y)S^u x (yMy); 



lim 

x—>y 

e^ X v p R(x) UlM (x)R{x)R(y)<P(y)S^T> x U(yy, 
e^?\R(x)SP Ufl (x)R(x)R(y)<f>(y)DxK(y); 
e^ x v p Ii(x)u^ x S^R(xMy)^yMy); 

e^\K(x)u»(x)D x K(x)ticf>(y)SWll{y) 



(56) 



By renormalization of loop-graph integrals, one gets a constraint on a linear combination of the UV-divergent parts 
of LEC's of the seven of the nine 2ir - 2nucleon contact counter terms of (^5|). First one calculates the amplitudes 
T corresponding to the contact terms of d55|). Then, writing on = a\ + XiL, (where L was defined in (|36])), and 
comparing with the UV-divergent portion of the total 1-loop amplitude ( ff<^ ) and Table II, one gets the following 
conditions (suppressing the two pairs of Pauli spinors and using ( [Tq ) - (|lj 

(1) spin-independent renormalization (graphs [(a)+...+(f)]): 



(V2) 2 ^4Ax + 2A 2 - 2A 3 J = (>/2) 
(2) spin-dependent renormalization (graphs [(g)+(h)]): 

(V2) 2 ^3A 5 + ^A 6 + 6A 7 -2A 9 



which on addition gives: 



(V2) 2 ( 4Ai + 2A 2 - 2A 3 + 3A 5 + -A 6 + 6A 7 - 2A< 



?A_JL 

12 128 



-(V2) 6 f , 



(V2)V15, 45 , 



24 V 16 



(57) 



(58) 



(59) 
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Note that the A4, As terms do not contribute at the kinematic point (|l7|). So, it is clear that the renormalization of 
the loop integrals can be accomplished by the 27r-2nucleon contact terms. 

Because the nucleon-nucleon interaction has a strong short range repulsion, the nucleon-nucleon wave function 
vanishes at short relative distances. Since the contact terms in and (^) behave as S^ 3S, (x — y), one does not get 
any contribution from these 2tt - 2nucleon contact terms, as well as from the UV-divergent parts of the 1-loop integrals. 
Hence, only UV-fmite parts of the loop integrals contribute to the DCX loop amplitudes after renormalization. 



VI. SUMMARY AND DISCUSSION 



The goal of this paper has been to calculate the one loop correction to the 2-nucleon amplitudes for pion double 
charge exchange (DCX) scattering by a nuclear target, at threshold, in the framework of HBChPT. For a numerical 
estimate of the 2-nucleon- 1-loop correction, the amplitudes are evaluated at a typical kinematic point ([l?]). These are 
compared to the leading order amplitudes (for the tree graphs), obtained in the same theoretical framework. 

The 2-nucleon- 1-loop graphs give a threshold contribution of about 0.016 relative to the leading order tree graphs at 
the same kinematics [with (p'—p) 2 = M 2 }. Both calculations are done in the static limit and the impulse approximation 
for the nucleons, and using the renormalized values of the axial-vector coupling, pion-decay constants and the pion 
and nucleon masses. As remarked at the end of Section V, these corrections are of the order of M 2 / '(A-kF^) 2 ~ 0.014 
times the tree graphs in LO (y — 0). This is as expected for a chiral correction of v = 2. 

As noted in Section II, there are also corrections to the LO tree graphs to order v = 2 arising from corrections to 
the 7r — 7r and 7r — N vertices, which have been studied in the literature. The 7r — 7r vertex has been determined to 
the required order (to one loop) 0], with the LECs fixed phenomenologically. However, although the n — N vertex 
corrections have also been studied through one loop 0,0, the LECs have only been determined for the on-shell 
threshold cases of 7r — N scattering lengths Q] and n^ 1 production on a single nucleon Jll]]. Therefore, there is not 
enough information from these studies to fix all the vertex LECs required to carry out the full v = 2 calculation of 
the DCX tree graphs. 

However, the 1-loop calculations cited do show large corrections to the vertices. We can estimate how much this 
would change the DCX tree-graph amplitudes by using the on-shell results as follows. First, we characterize the three 
amplitudes of (Bl]),(E3j) and fl29) in terms of vertex coefficients: 



T[ a ] ~ DxgA, 

T [b] ~ Ag% (60) 
T [c] ~ (a-) 2 . 

The amplitude D\ is the irN — > irirN threshold amplitude |ll| that contributes to the DCX tree graph [a]. To 
leading order, Di(LO) ~ gA/F 2 . The corrections through one loop give JllJ 

D x = £> 1 (LO)(1 + 0.67). (61) 

The loop correction is actually 0.15; the larger correction, 0.52, is from a 'recoil' term of order M 2 /m. 

The amplitude A gives the DCX contribution to tree graph [b] from the 7r — 7r vertex, and can easily be shown to 
be to the two s-wave scattering lengths a/, with isospins / = 0, 2, by 

A= (2a + a 2 )/3. (62) 

To LO, A(LO) ~ Mn/F 2 . The 1-loop corrections || change this to 

A = A(LO) (1 + 0.30). (63) 

The isovector ttN scattering length a - contributes to the double-scattering tree graph [c]. To LO, a~(LO) ~ 
Mtt/F 2 . With corrections through 1-loop, this becomes Q 

a~ = a~ (LO) (1 + 0.046). (64) 

This small correction masks a 20% contribution of the loop, partially cancelled by the 'recoil' correction. 

If we now use these corrected values for the vertices in the expressions (|60|), we obtain the following estimates of 
their effect to order v — 2 on the the tree-graph amplitudes: 
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T [a] = 1.67 T [a] (LO) 

T [b] =1.30 T [b] (LO) (65) 
T [c] = 1.092 T [c] {LO). 

These corrections are much larger than the 1-loop corrections of this paper, which involve both nucleons participating 
in the DCX reaction. However, the vertex corrections are only estimates, since not all the off-shell vertex coefficients 
are determined. 

So we are led to two conclusions about corrections to DCX to order v = 2. First, the one-loop correction to the 
tree graph amplitudes is small, whether compared to the LO tree amplitudes, or to those amplitudes also corrected 
to one-loop order. The dominant contributions (graphs 3b, 3c, and 4h) all have one nucleon propagator, and are 
proportional to g\ > 1. The largest, 4h, is clearly a correction to the largest tree graph, 2c, for double scattering. 
(Triple scattering on two nucleons has not been included; it belongs with other triple scattering graphs on three 
nucleons.) The result is weakly dependent on the choice of the renormalization point \i. 

For the purpose of renormalization of the loop graphs, we found that nine 27r - 2nucleon contact terms are required, 
of which seven contribute at the specific kinematic point considered (|l7|). However, since we have assumed that the 
nuclear wavefunctions vanish near zero (relative) distance of nucleon separation, none of the (UV-finite parts) of 
the contact terms will contribute, because of the (5' 3 '(x — y) in the forms of contact terms of ( pq ) and (|5^). If this 
anticorrelation assumption were relaxed, then one would have nine unknown LECs for ( |55| ) (seven for the kinematic 
point (p7|)), and five for (|5^), to be determined empirically, to the order considered [0(<r)J. 

The second conclusion is that the one-loop corrections to the tt — tt and tt—N vertices produce the largest corrections 
at order v = 2 to the LO tree amplitudes. The large correction to the tt — tt vertex has been related to the pion radius 
that is, to the correction for the LO assumption of a point pion. Presumably, similar effects, including resonant 
contributions, give the large tt — N vertex corrections. It is clear that these corrections are anomalously large, while 
the one-loop corrections to the two-nucleon amplitudes are of expected size for the 0(q 2 ) terms. The latter may reflect 
the fact that the N — N system is already of finite size. 

However, as noted above, the vertex corrections for the tt — N amplitudes have not been fixed sufficiently for a 
complete calculation of the DCX tree amplitudes. This would require an extension of the 1-loop theory of the various 
tt — N amplitudes to off-shell kinematics, beyond what has already been done in the literature. As a practical matter, 
the tt — N on-shell constants could be taken directly from experiment, as in the impulse approximation, rather than 
from higher orders in chiral perturbation theory. 
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APPENDIX 

APPENDIX A: VERTICES AND F'S FOR 1-LOOP GRAPHS FOR tt-NN 

The pion field is represented as a matrix-valued SU(2) field U defined as following: 

e 1 ^ = U(tt), (Al) 
where </> is a traceless (as det U=l) matrix written in terms of the pion triplet as: 

The matrix is written in a space that corresponds to the standard (p/n) isospinors which are implied in the H notation. 
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1. Vertices 



The leading order terms in (H)ChPT, written using (Al) and (A2), that are used for the evaluation of mir — NN(m = 
1, 2, 3, 4) and 4-7r vertices, are written out. 

IttNN : g A HS • u (1) H 



-(V2)^ 
fix 



-^pS • dn°p + pS • <97r + n + nS • dir'n - -^nS • dir°n 
v 2 v 2 



2ttNN : iRv ■ D (2) H = Hv ■ T (2) H 



PP7T 1 V 



+ V2pn7r [0 v • dir +] + \/2np7r [ "v • <9tt 01 



3ttNN : g A RS ■ u (3) H 
1 



= (V2) 



3 ffj 

12FJ 



71 1 



p 27r + 7r-S-9-7r u -7r u (7r + S-a7r-+7r-S-a7r+) p 



p( 7r° 2 S • dn+ - 7r+7r°S ■ Stt + TrV'S ■ Stt+I )n + h.c 



+ -^=n^-2vr + 7r-S • dir° + tt° + t^tt^S • dn + + tt+S • dn~) )n 
4ttNN : Eiv ■ D (4) H = Hiv • T (4) H 



(V2) 4 



96F£\ 2 
+nn7r^ + w ■ dir~^ 



+ 71" 7T 



pp7T [ ~v • dn +] + (v^pW+v • dn 0] + h.c.) 



F 2 / \ (4) 

4tt: ^-l(d^d fl U)+M*((U* + U-2)) 



= (V2)< 



1 



24FJ 



-27r+7r-^7r+^7r- + (^) 2 (a Al7 r+) 2 + (t^) 2 ^-) 2 + M 2 (7r+) 2 (7r-) ; 



(A3) 



2. Combinatoric Factors for 1-Loop Graphs 

The combinatoric factors "f of (Il9| ) for 1-loop graphs, are evaluated. First, (See (|l9|) ) is written out in terms 
of "H</> m H," m = 1,2,3,4 and "t^'v which represent the run — Uj nucleon vertices with m = 1,2,3, 4; = 2 and 
m = 4; rij = 0. Then, using the time ordering properties of bosonic and fermionic bilinear fields, those terms relevant 
to the DCX 1-loop graphs are picked out. Then on comparison with (|l9|), the combinatoric factors are read off for 
the eight 2-nucleon 1-loop graphs. They are listed along with the corresponding forms of equation (|n]) in Table V. 
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TABLE V. Combinatoric Factors for the 1-Loop Graphs 



s£/ = S/IIL 1 <* 4 ** 

^(nr =1 E,4(^)) 



^(nr =1 4(^)) 

relevant to x 



ir 



^lloil) ;x : .f llo 'll ) 



H(^ + ^ 2 + 3 )H (xi) 



H</> 3 H (x 3 



nti( fl (0+^)H 



+</) 4 (a* 



^H0hJ (xi)^H0H )(x 2 ) 



H0 2 H (x 3 )<?!> 4 (x4) 



4 nu^r 



H0 2 H (xi 



t~ 1 1, .f'i'r'r 



H^H (xi 



H<^. 2 H (x 2 ) H0 2 H (x 3 



H0 2 H (x 2 ) H0 2 H (x 3 ) 



H0 2 H + 4 (x 



II" 



H^H (xi) 



H0 2 H (x 2 )0 4 (x 3 ) 



H(0 2 + 4 )H (xi 



H0 2 H (xi) 



H(0 2 + 4 )H (x 2 ) 



H0 4 H (x 2 ) 



2 



nLr/d 4 ^ H0 3 H ( Xl 



H<£ 3 H (x 2 ) 



£ tt2 

2 



H0 3 H (xi) 



H<£ 3 H (x 2 ) 



n- =1 



; / d 4 ^ / d 4 x 2 / d 4 x 3 T 



fi^H j (xi) 



H(0 + <^ 2 + 3 )H ( Xi 



Ll.rll j:x,:^llo : ll 
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APPENDIX B: 1-LOOP INTEGRALS 



There are 11 1-loop integrals defined and evaluated in this appendix. The integrals 74^ do not occur in the 1-loop 
DCX amplitudes as discussed later. The notations used, though similar to the ones used by || and |fl3|| , are slightly 
different. 

Four of those integrals are referred to as basic integrals : A^O, M 2 ), J wN , J"* 71 , 70, evaluated at different kinematic 
points. (The in A^O, M 2 ) implies that the "external" 4-momentum squared "p 2 " that appears in other integrals, 
is zero) They are defined below. 



1 f d d k 1 , , 

ij WT M 2 -k 2 -^ A ^ M) 
1 r d d k 1 

~i J (2ir) d {vk-uj- ie)(M 2 — k 2 — ie) 
1 f d d k 1 



- J wn (V 2 ;M 2 ) 



i J (27r) d ) (M 2 - fc 2 - ie)(M 2 - (k - V) 2 - ie) 

1 r d d k 1 

1 J (2ir) d (vk-uj- ie)(M 2 - k 2 - ie)(M 2 - (fc - V) 2 - ie) 

= 70(^,0, V 2 ), (Bl) 

with fi = v ■ V. In [Of, V 2 that figures in J^^iV 2 , M 2 ) is > 0. In this paper, one requires V 2 < 0. 

The integrals A^M 2 ), J™{J> 2 \ M 2 ) and J wN (u) have already been evaluated in the literature ( HQ), where 
they were denoted by A ff (0), J(s or t) and Jq(cj), respectively. We have alterred the notation for A T and J to 
introduce M 2 as a variable, because we need to differentiate these expressions with respect to M 2 , as in (Bll) and 
(B20). The form of 7o(w, fi, V 2 ), required for DCX, has not been evaluated previously. It is a generalization of the 
integral denoted by jo(u>) in 0, for which tu = 0,V 2 = 0, V — > —V (SI used in this paper is denoted by "w" in p]). 
The reason for introducing 17 and P 2 in addition to u> in the argument of 70 becomes clear from the points (a) and 
(b) in the paragraph after (B3). Also, note that the indexing (i = 1,2,..) of J™(V 2 ] M 2 ) and ^iUu,Sl,V 2 ) in (Jb|) 
and (B3) below differs from that used for the corresponding integrals J%?(s or t) and 7»(a>) in |^Jl^|. 



The above four integrals are referred to as "basic" because the integrands have no momentum dependence in their 
numerators, and are the most basic of In propagator, l-nucleon-l7r propagators, 27r propagators and l-nucleon-27r 
propagators integrals respectively. 

The remaining integrals have momentum dependence in numerators (and hence have a tensorial character), and 
are defined as: 



1 f d k k^i or k^i ku 



7WCP 2 ;M 2 



or 



(27r) d ) (M 2 - fc 2 - ie)(M 2 - (fc - V) 2 - ie) 

g^JFP 2 ; m2 ) + *PJ>uJFP 2 \ M ') ( B2 ) 



1 f drk fc M or k^k^ 

1 J (2ir) d (vk-u- ie)(M 2 - fc 2 - ie)(M 2 - (fc - Vf - ie) 

^v^u^V^+V^u^V 2 ); 

or 5^73^, Si, V 2 ) + v^ 4 (lo, SI, V 2 ) + 7yVy&(w, 0, V 2 ) 

+ {^+^'PM^M,V 2 ), (B3) 

where alternative numerators in fc M or k^k v have been shown in each case. 

Other than notational differences, there are also the following kinematical differences between the the integrals of 
(B3) and those figuring in pj : (a) Our V 2 does not vanish, in general, unlike the case for fj], in which V 2 is the 
four-momentum squared of an on-shell photon; (h)u> ^ 0; SI = v ■ V ^ uj; (c) we consider (fc — V) 2 instead of (fc + V) 2 
in the integrals. 

Of the nine, two do not contribute to DCX: 74^. The reason is that because these two integrals can occur only in 
graphs (g) and (h) with integrands having numerators of the type • fcS' 2 ' • fc, their coefficients will have at least 
one S^ 1 )''- 2 -' • v, which vanishes in the static limit. (That is because in the static limit v 1 * = (1, 0) and = (0, ?)) The 
remaining seven can be rewritten as linear combinations of the four basic integrals of (Bl). 
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All the integrals for the calculation are real, and hence so are all the 1-loop amplitudes. This is because the elastic 
scattering amplitude becomes real at threshold. 

Two further simplifications of the calculation are obtained by the following. 
(A) 



The 37r-propagator integrals that occur in 1-loop graph Fig 7.2e using (16), are of the type 



d d k 



which can rewritten as 



i J (2ir) d (Af 2 -k 2 - ie) 2 (M% - (k - V) 2 - ie) ' 

d i r d d k /(/y 



dM 2 i J {2n) d (M 2 - k 2 - ie)(M% - (k - V) 2 - ie) ' 



(B4) 



(B5) 



where M 2 is set equal to M 2 after differentiation. 
(B) 

Use is also made of some identities involving the loop integrals, which are valid only in the static limit and impulse 
approximation (|l6|): 

7o (w = ±M v ,n = ±AU,V 2 = ( Pi -p 2 ±q^ 2 ) 2 ) 
= ^o(u J = 0,n = TM 7T ,V 2 = ( P 2-P4Tqi,2) 2 ) 
= lo (iu = 0,n = TM 7r ,V 2 = ( Pi ~p 2T q h2 ) 2 ) 
= 7o (u; = ±M 7r ,rj = ±A^,P 2 = (p 2 -p 4 ±q h2 ) 2 ), 



where p% « (0, -p), p\ « (0, -&), v» = (1, 0), v ■ S 1 ' 2 = 0. 
Also we uses the following the identity: 



7o (w = ±M W , n = ±M W , V 2 = M 2 - P 2 ) 
= 7o (w = 0,0 = TM n ,V 2 = M 2 - P 2 ). 



(B6) 



(B7) 



1. Basic Integrals 



The results after evaluation of the basic integrals of (Bl) are given below. 

M 2 M 2 

A„(0, M 2 ) = 2M 2 L + -—In— + 0(d - 4); 



167T 2 (I 2 

f2 

2 



V 2 )J™{V 2 ;M 2 ) = -2L 



p* Aw 2 
1 



.c'-nrcrosl | + £>(d-4); 



, M 2 (A-V 2 ), M 2 
1 + In^r - - Tr^ln- 



A?(V 2 ;M 2 ,M 2 ) 
-2V 2 



(Ai - V ) 



2^2 _ A 2 1 



(A3 + V 2 ) 2 - A 2 ) 



(A = M - M 2 - A(V 2 ;M 2 , M 2 ) = [V - (M - M^) 2 ][P Z - (M + M„) 2 ]) 



h=J™(~M 2 ,M 2 



-2L 



16tt 2 



M VV5-1 



(B8) 
(B9) 



I 8 = J^{V 2 = 0; M = M 2 ) = -2L 



167T 



M 2 

(l + ln -) + 0(d-4) 



(B10) 
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0( ~ p2) ^^ (p2;M2 ' M * )IM2 = M * 

1/1 1 



16^ 2 \M% v /(-P 2 )(4M2-p2) 



In 



^4Af2 - -P 2 + ^/T^2 
v/4M2 - P2 _ ^/T^2 



(*) 7 o(^ = 0, fi = M T ,P 2 = 0) - 7o (a; = 0, ft = -M„ P 2 = 0) = - 

(m) J n = 7o (w = -M ff ,fi = 0,P 2 = — M 2 ) 
Zn 



32M W ' 



(BH) 
(B12) 



k 1 + \ f^) -4^ \r-^) 

i , ri + v^r 



SttM^ 



-Zri 



(B13) 



in is a new integral that one needs to evaluate for the DCX 1-loop graphs. For values of V 2 + 1 0, — M 2 , one gets the 
following expressions, which are also new: 

(i)70(w = -M ff) ft = -M^P 2 = M 2 - P 2 ) 
-7o(w = M„ n = M ff , P 2 - M 2 - P 2 ) 



where 



1 



2n\P\ 
\P\ 



— arctan 



E 



M 2 +P 2 

(2k)\ 



M 77^^ 2-(fc!) 2 (2fc + l) 

/3 P 2 M 2 

xFi -,fc+ 1,1,2;-— 7 , 

\2' M2'[2M2 + P2] 



Fi(a, 6, c, d; y) = Appell functions(i — 1, 2, 3, 4); 



Fi(a,6,c,d;a;,2/) = E E 



m— n— 



c m+ „ m!n! 



(B14) 



(B15) 



(i*)7o(w = -m W) ft = o, v 2 = -p 2 ) = 

(2k)\ 



1 + 



p 2 _ ]P| 

2M? Mtt 



P 2 



16tt 2 |P|/ 



4 + 



P 2 



E 

fc=0 



{k\) 2 2 2k (2k + l) 



25^ + 2,^2^^ + 2,1;' 



p^ 

_\W 2 



4 + 



P 2 



-a+^)+vM| 



-B(fc + l,l) 2^^ + 1,1;^+ |; 





^ 2 + /a 4. -f 32 




-a+^+^+n 
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167r 2 |PU/4 



p 2 

M 2 

- 1 , x k (2x- 1) 
< / — -=-- 
Jo Vl-x 



32 



fe=0 



(2fc)! 
(fc!) 2 2 2fc (2A: + 1) 



j32 P 



(B16) 



2. Other Integrals 

The results obtained after evaluation of other integrals of (52) — (-B3) are given below, 
(a) 2tt — Propagator Integrals 

The following relations, (B20), (B22), (B24), (B26) and (B28), are not explicitly given in the literature. These 
relations are absolutely general and become relevant when evaluating the off-threshold 1-loop amplitudes for kinematic 
points away from fll^). 



i 9 = jr(o,M^) = --A 7I (M',oy 1 



(B17) 



ho = JT%M%) = \j™Q,M$) = \h 



(B18) 



J™{V 2 ,M Z ) 



1 

2^2 



I 2 = jr{-MlMl) = -J™(-M 2 ,M 2 ) 



A T (0, M 2 ) - A w (0, Ml) + (A + T 2 )J{T 2 ,M 2 ) 



1 



32tt 2 



-1 + In- 



Ml 



{ V5 + 1 
\V5-1 



(B19) 



2 ^ 

-M„- T ^ T Ji 7r \ M 2 =M 2. V 2 = _ M 2 



d 



1 



2 V dM 2 



A T (0,Af 2 )-Af2 



<9M 2 



1 



J™(-M 2 ,M 2 ) 



\/5 , /^V5 + 1 



2 32tt 2 40tt 2 Va/5-1 



M 2 =M2:T> 2 =-M 2 



(B20) 

(B21) 



2M. 



f2 



J^{V 2 ;M 2 =M 2 ) = 

which for V 2 = -M 2 gives : 
h = J^{-Ml,Ml) 



Ei 

2 



j™^) - A 7r (0,M 2 



2(d-l) 



6 



2!) 



576tH 



-Mi 



192tH 



rMiZn- 



5-y/5 , ^V5 + 1 
V5- 1 



192?H 



An 



O(d-A); 



(B22) 



(B23) 



24 



J™(T 2 ;M 2 = M*) = J^(P ) 

dV 2 - AM 2 J J^iV 2 ) + 2(2 - d)A*(0, M|) 



4(d- 1)P 2 
which for P 2 = — M 2 gives : 

i 4 ^jr(-M 2 ,M 2 ) 
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= 4i 



1 . Ml 
rln- 



?(-p 2 ^ 



dM 2 
1 



3 288tt 2 48tt 2 /i 2 24tt 2 

0(d-4). 
J^{V 2 ;M 2 ) 



2(1 - d) 

L 1 
2 ~ 192tt 2 



1 1 

. f4M 2 - V 2 \ — — T 7 ™ 



i a 



2<9M 2 



A ff (0,M 2 



r(l 



V 



M 2 ' 64tt 2 ^i 2 



1 /[4Af 2 -7> 2 ] 



64tt 2 
So, 



_p2 



in 



^4M 2 - V 2 + V^P 2 " 



^4Af 2 - 7> 2 - V^P 2 
i 1 



<9A/ 2 



jr(-M 2 ,M 2 )\ M 2 =M2 = - 



d 



0(d-4); 



1 , M 2 

An 



2 96tt 2 64tt 2 ^ 2 64tt 2 \ \/5 - 1 / 



So, 



/ 7 = 



<9Af 2 



9M 2 



J?"(-MZ,M£\ M 2 =M 2 =- 



1 =^4- 



96tt 2 16a/5 V-v/5-1 



(B24) 



(B25) 



(B26) 



(B27) 



(B28) 



(B29) 



(b) 2tt — INucleon Propagator Integrals 

The following relations, (B30), (B32) and (B35), are not explicitly given in the literature, and are absolutely general 
and become relevant when evaluating the off-threshold 1-loop amplitudes for kinematic points away from (|l7|). 



7l (^,ft,7> 2 ) = 



ZP*j**(p3) + -p 2 {2uj - fi)7o(o;, ft,P 2 ) -ft J* N (w) - J tA (w - 0) 



2CP 2 - ft 2 ) 

^(w) - - 0) - 2QJ™(V 2 ) + {V 2 - 2wft) 7o (w, ft, 7> 2 ) 



Hence, 



2(P 2 - ft 2 ) 

7i2 = 72(0, —M v , 0) - 72 (0, M ff , 0) = 



4M w Tr 2 



(B30) 



(B31) 



73 (^,ft,7> 2 ) = 



2M 2 7o (w, ft, "P 2 ) - 2cj 7i (w, ft, V 2 ) - P 2 72 (w, ft, V 2 ) - J* N (w - ft) 

2(d- 2) 
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Hence, 



Thus: 



Now, 



Hence, 



Thus: 



1 
4 

(2lu - Q) 



2M 2 7o (w, 0,P 2 ) - 2w 7l (w, n, P 2 ) - P 2 j 2 (uj, 0, P 2 ) - J Tjv (u; - 0) 
+ 0(d-4). 



32tt 2 



73 (w = 0, = -M^V 2 = M 2 - P 2 ) - 73 (co = 0,0 = M n , V 2 = M 2 - P 2 ) 



7 o(0, -M w , M 2 - P 2 ) - 70 (0, M-n, M 2 - P 2 ) 



+ [M y ] (j* N (-M„)-J" N (M„) 



[M 2 + P 2 
2P 2 



/ 13 = 73 (0, -M w , 0) - 73 (0, M„ 0) 
1 



7»(w,n, p 2 ) 



(d- 2)0,7™ (P 2 ) 



2(d-2)(P 2 -0 2 ) 
+2 72 (a), 0, P 2 )[(d - 1)P 2 + wQ(d - 2)] 

+(d - 3) J^w - fi) - 2w7i(w, 0, P 2 ) + 2M 2 7o (w, fi, P 2 ) 



2w7i(w,0,P 2 ) - + (3P 2 - 4wO) 72 (w,0,P 2 ) - J^w - 0) 

+ 0(d-4). 



4(P 2 - O 2 ) 



-2M 2 70 (-,O,P 2 )+ (2a; " } 



8tt 2 



7 5 (w = 0, = -M„ P 2 = M 2 - P 2 ) - 75(a) = 0, SI = M n , V 2 = M 2 - P 2 ) 



1 

4P 2 



[3M 2 - 5P 2 ] / J7lN 



+ 



2P 2 

(3[M 2 -P 2 ] 2 +4M 2 P 2 

~ 2P 2 
2M 7r [5M 2 -3P 21 



7 o(0, -M x , M 2 - P 2 ) - 7o (0, M,r, M 2 - P 2 ) 



P 2 



^(M 2 -P2 ) + ^| 



0(d-4). 



J14 = 75(0, -M ff , 0) - 75(0, Mtt, 0) 
= l^ + 64k + ° ( ^ 4) - 



(B32) 



(B33) 



(B34) 



(B35) 



(B36) 



(B37) 



2G 



[1] S.R.Beane, C.Y.Lee and U.van Kolck, Phys. Rev. C 52, 2914 (1995). 

[2] B-Y Park, F.Myhrer, J.R.Morones, T.Meissner and K.Kubodera, Phys. Rev. C 53, 1519 (1996). 

[3] J.Gasser and H.Leutwyler, Annals of Physics 158, 142 (1984). 

[4] V. Bernard, N. Kaiser and Ulf-G Meissner, Phys. Lett. B 309, 421 (1993). 

[5] M.F.Jiang and D.S.Koltun, Phys. Rev.C 42, 2662 (1990). 

[6] V.Bernard, N.Kaiser and Ulf-G.Meissner, Int. J. Mod. Phys. E4, 193 (1995). 

[7] S.Weinberg, Phys. Lett. B 251, 288 (1990). 

[8] S.Weinberg, Phys. Lett. B295, 114 (1992). 

[9] David B. Kaplan, Martin J. Savage, Mark B. Wise, Phys. Lett. B 424, 390 (1998). 
[10] David B. Kaplan, Martin J. Savage, Mark B. Wise, Nucl. Phys. B 534, 329 (1998). 
[11] V.Bernard, N.Kaiser and Ulf-G Meissner, Nucl.Phys.B 457, 147 (1995). 
[12] M.G.Olsson and L.Turner, Phys. Rev. Lett 20, 1127 (1968). 
[13] J.Gasser, M.E.Sainio and A.Svarc, Nucl. Phys.B 307, 779 (1988). 



(i) 



(ii) 



(iii) 



(iv) 



(v) 

FIG. 1. Elementary Vertices 



27 



(a) 



(b) 



(c) 



FIG. 2. Tree Graphs 



(a)(1) 



= 



(a)(2) 



(b)(1) 



(b)(2) 



(c)(1) 



(c)(2) 



(d)(1) 



(d)(2) 



FIG. 3. 2Nucleon-l Loop Graphs (a) -(d) 
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(e) 



(f) 



(g)(1) 



(g)(2) 



(h)(1) 



(h)(2) 



FIG. 4. 2Nucleon-l Loop Graphs (e) - (h) 
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FIG. 5. 2tt-2 Nucleoli Contact Graph 
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